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pL^ ' Abstract. Close-in extrasolar giant planets may be directly detectable 

\^ , by their reflected light, due to the proximity of the planet to the illu- 

fN| ' minating star. The spectrum of the system will contain a reflected light 

component that varies in amplitude and Doppler shift as the planet orbits 
the star. Intensive searches for this effect have been carried out for only 
one extrasolar planet system, r Boo. There exist several other attractive 
qq ' targets, including the transiting planet system HD 209458. 
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Introduction and Motivation 



Radial velocity surveys of nearby Sun-like stars have uncovered a population of 
close-in orbiting companions of roughly Jupiter mass. The ten such objects with 
semi-major axes a <^ 0.1 AU are listed in Table 1, along with the values for the 
semi-major axes, (minimum) masses, spectral types of the stars, and equilibrium 
temperatures (calculated from the estimated values for the semi-major axes, 
stellar radii and effective temperatures, and assuming a Bond albedo of ^4). 

5_i ■ Table 1. Close-in Planets Ordered by Increasing Semi-Major Axis 



Star 


Spectral Type 


a (AU) 


M p (M Jup ) 


T cq (l - A)~ 


HD 187123 


G3 V 


0.042 


0.52/sini 


1400 


HD 75289 


GO V 


0.046 


0.42/sini 


1600 


r Boo 


F7 V 


0.046 


3.87/sini 


1600 


HD 209458 


GO V 


0.047 


0.69 


1500 


51 Peg 


G2 V 


0.051 


0.45/sin i 


1300 


v And 


F8 V 


0.059 


0.71/sini 


1500 


HD 217107 


G7 V 


0.072 


1.28/sini 


1000 


HD 130322 


K0 V 


0.088 


1.08/sini 


800 


55 Cnc 


G8 V 


0.11 


0.84/sini 


700 


Gl 86 


K0 V 


0.11 


3.6/sin i 


700 



A successful spectroscopic detection of an extrasolar planet in reflected light 
would yield the inclination, and hence the planetary mass, and would also mea- 
sure a combination of the planetary radius and albedo. Furthermore, it would 
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open the way to direct investigation of the spectrum of the planet itself. Con- 
versely, a low enough upper limit would provide useful constraints on the radius 
and albedo of the companion. 

The predicted albedo of a close-in extrasolar giant planet has been the focus 
of recent theoretical work (Marley et al. 1999; Seager, Whitney, & Sasselov 
2000), and the possible values range by several orders of magnitude based on 
the atmospheric temperature, chemical composition (see, e.g., Burrows & Sharp 
1999), and the presence (or absence) of atmospheric condensates, and their 
respective size distributions. 

2. The Method of Reflected Light 

The amplitude, relative to the star, of the observed flux reflected from a close-in 
planet viewed at a phase angle a (where a is the angle between the star and the 
observer as viewed from the planet) is given by 



where R p is the planetary radius, a is the semi-major axis, and p\ is the 
wavelength-dependent geometric albedo. The quantity e is the contrast ratio 
at opposition. The phase function <p\{a) is the brightness of the planet viewed 
at an angle a relative to its value at opposition (where by definition (f)\(0) = 1), 
and is a monotonically decreasing function of a for all physically reasonable 
atmospheres. For planets detected by the radial velocity technique alone, a is 
determined, and a at a given time is prescribed by the known orbital phase $ 
and the unknown orbital inclination i, but R p , p\ and (f>\(a) remain unknown. 

A distant observer viewing the solar system in the V band would measure 
a tiny reflected light ratio of / ~ 4 x 10 -9 from Jupiter viewed at opposition. 
This ratio increases dramatically for close-in, giant planets, due to the small 
orbital separation. Scaling equation (1), we find 



This photometric modulation is most likely beyond the reach of ground-based 
observations, but will be accessible to upcoming photometric satellite missions, 
such as COROT (see, e.g., Rouan et al. (1999)), provided such missions can 
achieve a precision of ~ 20 /xmag with stability over timescales of a few days. 

Photometry alone will not yield the orbital inclination, and hence the plan- 
etary mass will remain constrained only by the lower limit imposed by the radial 
velocity observations. However, a spectrum of a close-in planet system will con- 
tain a secondary component that varies in brightness according to equation (1) 
and in Doppler shift (relative to the star) according to 




(1) 



f x (a) ~ 9.1 x 10- 5 Px 




(2) 



v p (<S>) = -K s 



M s + M p 



COS 2lT& 



(3) 



where M s and M p are the stellar and planetary masses. The stellar radial veloc- 
ity amplitude K s and orbital phase 3> are determined from the observed radial 
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Figure 1. The solid curves are the confidence levels on the upper 
limit for the relative reflected flux e in the r Boo system as a func- 
tion of orbital inclination i, if the reflected light is a copy of the stellar 
spectrum. The dashed curves are the same confidence levels under the 
assumption that the system is tidally locked and thus the planet reflects 
a non-rotationally-broadened copy of the stellar spectrum. Upper lim- 
its on the geometric albedo p under the assumption that R p = 1.2 Rj up 
are shown on the right-hand axis, and the values for Jupiter and Uranus 
are included for comparison. The lack of transits (Baliunas et al. 1997) 
excludes i ^ 83°, and i <J 17° can be excluded under the assumption 
that the stellar rotation axis is co-aligned with that of orbital motion, 
as it would imply an abnormally large rotational velocity (Gray 1982). 



velocity orbit of the star. For the close-in planets, \v p \ ~ 100 kms , which is 
two orders of magnitude larger than the resolution of current spectrographs. If 
the reflected light spectrum of the planet is detected, then the planetary mass 
is determined by equation (3). 

Furthermore, as first discussed in Charbonneau, Jha & Noyes (1998), if the 
star has been tidally spun-up so that its rotation period is equal to the plane- 
tary orbital period, then the planet would reflect a non-rotationally-broadened 
spectrum, further distinguishing these features from the stellar lines. 

In summary, the observational challenge is similar to that of transforming 
a single-lined binary system into a double-lined system, in the case of an excep- 
tionally large contrast ratio between the two components. See Charbonneau et 
al. (1999), hereafter C99, for more details. 
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3. Results for r Boo 

Figure 1 summarizes our search for the reflected light spectrum from the planet 
orbiting the star r Boo, as described in C99. At the 99% confidence level, we 
find no evidence for a reflected flux ratio in excess of 1 x 10 -4 . For edge-on 
values of the inclination (i ^ 70°), this ratio is further restricted to be less than 
5 x 1CT 5 . Assuming a planetary radius of 1.2 i?j up (Guillot et al. 1996), this 
limits the geometric albedo to p\ < 0.3 for 466 nm < A < 498 nm. 

These conclusions are marginally in conflict with the recent claim of a pos- 
sible detection of reflected light for the same system by Cameron et al. (1999). 
They find a reflected light ratio of (1.9 =L 0.4) x 10~ 4 in the wavelength range 
from 456 nm to 524 nm, and do not detect the planet in wavelength bands on 
either side of this region (385 nm to 456 nm, and 524 nm to 611 nm). Figure 2 
demonstrates the confidence intervals that we would have obtained on the two 
unknowns {e, i} for inclinations near the one claimed by Cameron et al. (1999) 
and for a selection of contrast ratios. 

Taking into account the different descriptions of the assumed phase function 
4>\(a), and the respective errors stated in the two papers, it is possible that the 
two results are consistent, i.e. the inclination claimed by Cameron et al. (1999) is 
correct, but the amplitude of the signal is less than their most probable value by 
1.5 times their quoted error. However, Cameron et al. (1999) caution that there 
is a 5% possibility that their detection is spurious. This controversy should be 
resolved by observations in the near future that can now be tailored (in selection 
of wavelength region and orbital phase) to verify the claimed detection. 

4. Additional Targets 

There are two principal considerations that enter into the choice of which planets 
may be most easily detected by reflected light: An ideal system is one with a 
small semi-major axis and a large stellar apparent brightness, and thus a large 
reflected light signal relative to the photon noise of the star. There are two 
other systems that satisfy these criteria, and we compare below the prospects 
for detecting each of these with those for r Boo (a = 0.046 AU, V = 4.5). 

The planet orbiting v And is further from its star (a = 0.059 AU), but the 
star has a greater apparent brightness (V = 4.1). Combining these two effects, 
the detection threshold that could be established (given an equivalent amount 
of observing time) would be 1.36 times higher relative to the primary star than 
that for r Boo. For the planet orbiting 51 Peg, the larger semi-major axis 
(a = 0.051 AU) and fainter star (V = 5.5) result in a detection threshold that is 
1.93 times higher than that for r Boo. 

Near edge-on inclinations {i ~ 90°) are desirable for reflected light observa- 
tions since they allow the planet to be viewed at the full range of phase angles. 
There are additional considerations that may constrain the inclination for these 
two systems. Assuming that the orbital planes of the three planets orbiting 
v And (Butler et al. 1999) are co-aligned, stability arguments (Laughlin & 
Adams 1999) favor suu ^ 0.75. However, Hipparcos astrometry favors (Mazeh 
et al. 1999) sini ~ 0.4. The rotation periods of v And and 51 Peg indicate that 
the stellar rotation has not become tidally locked to the orbital period, and this 
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Figure 2. The detection thresholds shown in Figure 1 were tested by 
directly injecting a secondary into the data at a given amplitude and 
inclination. Three such tests are shown above, and the black crosses 
indicate the input parameters. The contours demarcate the 68% and 
95% confidence regions on the derived parameters for the planet. 

in turn might imply planetary masses which are near the minimum mass values 
(Drake et al. 1998). 

We expect 51 Peg b and v And b to have larger radii than the more massive 
r Boo b (Guillot et al. 1996). If r Boo b has been detected, it has a very large 
surface area. If 51 Peg b and v And b have comparable albedos to r Boo b, they 
should have a comparable or larger reflected light ratio. 

5. The Transiting System HD 209458 

In the case of a planet that is observed to transit its star, the situation improves 
considerably. Since transit observations yield R p and i, then only remaining 
unknowns in equation (1) are p\ and <f>\ (a). 

The first transiting extrasolar planet, HD 209458 b, was detected by Char- 
bonneau et al. (2000) and Henry et al. (2000). The dominant uncertainty in 
deriving the planetary radius and orbital inclination is the estimation of the 
stellar radius, and not the photometric precision. A detailed analysis of the 
stellar modeling, and the resulting derivation of the planetary parameters and 
the systematic uncertainties is presented in Mazeh et al. (2000). They derive 
R p = 1.40 ± 0.17 i? Jup , a = 0.467 AU, and % = 86? 1 ± 1?6. Substituting these 
values into equation (3), we find fx (a) = 2.1 x 10~ 4 p\ 4>\{a). 
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Transiting systems such as HD 209458 are particularly desirable targets for 
reflected light measurements since by securing measurements near opposition 
(where <p\ (0) ~ 1), an observed reflected light amplitude is a direct measurement 
of the geometric albedo. Furthermore, one can use the spectra obtained while 
the planet is behind the star (rather than near inferior conjunction) to create 
the required stellar template spectrum. One difficulty is that HD 209458 is 
significantly fainter (V = 7.6) than the systems described above. 

We have simulated what could be achieved in a modest observing run using 
the Keck-1 telescope and HIRES spectrograph (Vogt et al. 1994), based on 
experience gained from our earlier observing run on r Boo (C99). In three nights 
when the planet is near opposition, we would obtain a detection threshold of 
4 x 10~ 5 for the reflected flux ratio. Thus we would constrain the geometric 
albedo to be p\ < 0.2, or directly detect the planet in reflected light. 
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